Reduced Protein Stability of 11 Pathogenic Missense STXBP1/MUNC18-1 Variants and Improved Disease Prediction.

BACKGROUND: Pathogenic variants in STXBP1/MUNC18-1 cause severe encephalopathies that are among the most common in genetic neurodevelopmental disorders. Different molecular disease mechanisms have been proposed, and pathogenicity prediction is limited. In this study, we aimed to define a generalized disease concept for STXBP1-related disorders and improve prediction. METHODS: A cohort of 11 disease-associated and 5 neutral variants (detected in healthy individuals) were tested in 3 cell-free assays and in heterologous cells and primary neurons. Protein aggregation was tested using gel filtration and Triton X-100 insolubility. PRESR (predicting STXBP1-related disorder), a machine learning algorithm that uses both sequence- and 3-dimensional structure-based features, was developed to improve pathogenicity prediction using 231 known disease-associated variants and comparison to our experimental data. RESULTS: Disease-associated variants, but none of the neutral variants, produced reduced protein levels. Cell-free assays demonstrated directly that disease-associated variants have reduced thermostability, with most variants denaturing around body temperature. In addition, most disease-associated variants impaired SNARE-mediated membrane fusion in a reconstituted assay. Aggregation/insolubility was observed for none of the variants in vitro or in neurons. PRESR outperformed existing tools substantially: Matthews correlation coefficient = 0.71 versus <0.55. CONCLUSIONS: These data establish intrinsic protein instability as the generalizable, primary cause for STXBP1-related disorders and show that protein-specific ortholog and 3-dimensional information improve disease prediction. PRESR is a publicly available diagnostic tool.

Multi-omic and functional analysis for classification and treatment of sarcomas with FUS-TFCP2 or EWSR1-TFCP2 fusions.

Linking clinical multi-omics with mechanistic studies may improve the understanding of rare cancers. We leverage two precision oncology programs to investigate rhabdomyosarcoma with FUS/EWSR1-TFCP2 fusions, an orphan malignancy without effective therapies. All tumors exhibit outlier ALK expression, partly accompanied by intragenic deletions and aberrant splicing resulting in ALK variants that are oncogenic and sensitive to ALK inhibitors. Additionally, recurrent CKDN2A/MTAP co-deletions provide a rationale for PRMT5-targeted therapies. Functional studies show that FUS-TFCP2 blocks myogenic differentiation, induces transcription of ALK and truncated TERT, and inhibits DNA repair. Unlike other fusion-driven sarcomas, TFCP2-rearranged tumors exhibit genomic instability and signs of defective homologous recombination. DNA methylation profiling demonstrates a close relationship with undifferentiated sarcomas. In two patients, sarcoma was preceded by benign lesions carrying FUS-TFCP2, indicating stepwise sarcomagenesis. This study illustrates the potential of linking precision oncology with preclinical research to gain insight into the classification, pathogenesis, and therapeutic vulnerabilities of rare cancers.

DLG1 functions upstream of SDCCAG3 and IFT20 to control ciliary targeting of polycystin-2.

Polarized vesicular trafficking directs specific receptors and ion channels to cilia, but the underlying mechanisms are poorly understood. Here we describe a role for DLG1, a core component of the Scribble polarity complex, in regulating ciliary protein trafficking in kidney epithelial cells. Conditional knockout of Dlg1 in mouse kidney caused ciliary elongation and cystogenesis, and cell-based proximity labelling proteomics and fluorescence microscopy showed alterations in the ciliary proteome upon loss of DLG1. Specifically, the retromer-associated protein SDCCAG3, IFT20 and polycystin-2 (PC2) were reduced in cilia of DLG1 deficient cells compared to control cells. This phenotype was recapitulated in vivo and rescuable by re-expression of wildtype DLG1, but not a Congenital Anomalies of the Kidney and Urinary Tract (CAKUT)-associated DLG1 variant, p.T489R. Finally, biochemical approaches and Alpha Fold modelling suggested that SDCCAG3 and IFT20 form a complex that associates, at least indirectly, with DLG1. Our work identifies a key role for DLG1 in regulating ciliary protein composition and suggests that ciliary dysfunction of the p.T489R DLG1 variant may contribute to CAKUT.

De-Suppression of Mesenchymal Cell Identities and Variable Phenotypic Outcomes Associated with Knockout of Bbs1.

Bardet-Biedl syndrome (BBS) is an archetypal ciliopathy caused by dysfunction of primary cilia. BBS affects multiple tissues, including the kidney, eye and hypothalamic satiety response. Understanding pan-tissue mechanisms of pathogenesis versus those which are tissue-specific, as well as gauging their associated inter-individual variation owing to genetic background and stochastic processes, is of paramount importance in syndromology. The BBSome is a membrane-trafficking and intraflagellar transport (IFT) adaptor protein complex formed by eight BBS proteins, including BBS1, which is the most commonly mutated gene in BBS. To investigate disease pathogenesis, we generated a series of clonal renal collecting duct IMCD3 cell lines carrying defined biallelic nonsense or frameshift mutations in Bbs1, as well as a panel of matching wild-type CRISPR control clones. Using a phenotypic screen and an unbiased multi-omics approach, we note significant clonal variability for all assays, emphasising the importance of analysing panels of genetically defined clones. Our results suggest that BBS1 is required for the suppression of mesenchymal cell identities as the IMCD3 cell passage number increases. This was associated with a failure to express epithelial cell markers and tight junction formation, which was variable amongst clones. Transcriptomic analysis of hypothalamic preparations from BBS mutant mice, as well as BBS patient fibroblasts, suggested that dysregulation of epithelial-to-mesenchymal transition (EMT) genes is a general predisposing feature of BBS across tissues. Collectively, this work suggests that the dynamic stability of the BBSome is essential for the suppression of mesenchymal cell identities as epithelial cells differentiate.

Shifts in mutation spectra enhance access to beneficial mutations.

Biased mutation spectra are pervasive, with wide variation in the magnitude of mutational biases that influence genome evolution and adaptation. How do such diverse biases evolve? Our experiments show that changing the mutation spectrum allows populations to sample previously undersampled mutational space, including beneficial mutations. The resulting shift in the distribution of fitness effects is advantageous: Beneficial mutation supply and beneficial pleiotropy both increase, while deleterious load reduces. More broadly, simulations indicate that reducing or reversing the direction of a long-term bias is always selectively favored. Such changes in mutation bias can occur easily via altered function of DNA repair genes. A phylogenetic analysis shows that these genes are repeatedly gained and lost in bacterial lineages, leading to frequent bias shifts in opposite directions. Thus, shifts in mutation spectra may evolve under selection and can directly alter the outcome of adaptive evolution by facilitating access to beneficial mutations.

A targeted multi-proteomics approach generates a blueprint of the ciliary ubiquitinome.

Establishment and maintenance of the primary cilium as a signaling-competent organelle requires a high degree of fine tuning, which is at least in part achieved by a variety of post-translational modifications. One such modification is ubiquitination. The small and highly conserved ubiquitin protein possesses a unique versatility in regulating protein function via its ability to build mono and polyubiquitin chains onto target proteins. We aimed to take an unbiased approach to generate a comprehensive blueprint of the ciliary ubiquitinome by deploying a multi-proteomics approach using both ciliary-targeted ubiquitin affinity proteomics, as well as ubiquitin-binding domain-based proximity labelling in two different mammalian cell lines. This resulted in the identification of several key proteins involved in signaling, cytoskeletal remodeling and membrane and protein trafficking. Interestingly, using two different approaches in IMCD3 and RPE1 cells, respectively, we uncovered several novel mechanisms that regulate cilia function. In our IMCD3 proximity labeling cell line model, we found a highly enriched group of ESCRT-dependent clathrin-mediated endocytosis-related proteins, suggesting an important and novel role for this pathway in the regulation of ciliary homeostasis and function. In contrast, in RPE1 cells we found that several structural components of caveolae (CAV1, CAVIN1, and EHD2) were highly enriched in our cilia affinity proteomics screen. Consistently, the presence of caveolae at the ciliary pocket and ubiquitination of CAV1 specifically, were found likely to play a role in the regulation of ciliary length in these cells. Cilia length measurements demonstrated increased ciliary length in RPE1 cells stably expressing a ubiquitination impaired CAV1 mutant protein. Furthermore, live cell imaging in the same cells revealed decreased CAV1 protein turnover at the cilium as the possible cause for this phenotype. In conclusion, we have generated a comprehensive list of cilia-specific proteins that are subject to regulation via ubiquitination which can serve to further our understanding of cilia biology in health and disease.

Never-homozygous genetic variants in healthy populations are potential recessive disease candidates.

The rapid pace with which genetic variants are now being determined means there is a pressing need to understand how they affect biological systems. Variants from healthy individuals have previously been used to study blood groups or HLA diversity and to identify genes that can apparently be nonfunctional in healthy people. These studies and others have observed a lower than expected frequency of homozygous individuals for potentially deleterious alleles, which would suggest that several of these alleles can lead to recessive disorders. Here we exploited this principle to hunt for potential disease variants in genomes from healthy people. We identified at least 108 exclusively heterozygous variants with evidence for an impact on biological function. We discuss several examples of candidate variants/genes including CCDC8, PANK3, RHD and NLRP12. Overall, the results suggest there are many, comparatively frequent, potentially lethal or disease-causing variants lurking in healthy human populations.

Next Generation Protein Structure Predictions and Genetic Variant Interpretation.

The need to make sense of the thousands of genetic variants uncovered every day in terms of pathology or biological mechanism is acute. Many insights into how genetic changes impact protein function can be gleaned if three-dimensional structures of the associated proteins are available. The availability of a highly accurate method of predicting structures from amino acid sequences (e.g. Alphafold2) is thus potentially a great boost to those wanting to understand genetic changes. In this paper we discuss the current state of protein structures known for the human and other proteomes and how Alphafold2 might impact on variant interpretation efforts. For the human proteome in particular, the state of the available structural data suggests that the impact on variant interpretation might be less than anticipated. We also discuss additional efforts in structure prediction that could further aid the understanding of genetic variants.

Defining clinical subgroups and genotype-phenotype correlations in NBAS-associated disease across 110 patients.

PURPOSE: Pathogenic variants in neuroblastoma-amplified sequence (NBAS) cause an autosomal recessive disorder with a wide range of symptoms affecting liver, skeletal system, and brain, among others. There is a continuously growing number of patients but a lack of systematic and quantitative analysis. METHODS: Individuals with biallelic variants in NBAS were recruited within an international, multicenter study, including novel and previously published patients. Clinical variables were analyzed with log-linear models and visualized by mosaic plots; facial profiles were investigated via DeepGestalt. The structure of the NBAS protein was predicted using computational methods. RESULTS: One hundred ten individuals from 97 families with biallelic pathogenic NBAS variants were identified, including 26 novel patients with 19 previously unreported variants, giving a total number of 86 variants. Protein modeling redefined the ß-propeller domain of NBAS. Based on the localization of missense variants and in-frame deletions, three clinical subgroups arise that differ significantly regarding main clinical features and are directly related to the affected region of the NBAS protein: ß-propeller (combined phenotype), Sec39 (infantile liver failure syndrome type 2/ILFS2), and C-terminal (short stature, optic atrophy, and Pelger-Huët anomaly/SOPH). CONCLUSION: We define clinical subgroups of NBAS-associated disease that can guide patient management and point to domain-specific functions of NBAS.

Wobbling Forth and Drifting Back: The Evolutionary History and Impact of Bacterial tRNA Modifications.

Along with tRNAs, enzymes that modify anticodon bases are a key aspect of translation across the tree of life. tRNA modifications extend wobble pairing, allowing specific ("target") tRNAs to recognize multiple codons and cover for other ("nontarget") tRNAs, often improving translation efficiency and accuracy. However, the detailed evolutionary history and impact of tRNA modifying enzymes has not been analyzed. Using ancestral reconstruction of five tRNA modifications across 1093 bacteria, we show that most modifications were ancestral to eubacteria, but were repeatedly lost in many lineages. Most modification losses coincided with evolutionary shifts in nontarget tRNAs, often driven by increased bias in genomic GC and associated codon use, or by genome reduction. In turn, the loss of tRNA modifications stabilized otherwise highly dynamic tRNA gene repertoires. Our work thus traces the complex history of bacterial tRNA modifications, providing the first clear evidence for their role in the evolution of bacterial translation.

The Frequency of Internal Shine-Dalgarno-like Motifs in Prokaryotes.

In prokaryotes, translation initiation typically depends on complementary binding between a G-rich Shine-Dalgarno (SD) motif in the 5' untranslated region of mRNAs, and the 3' tail of the 16S ribosomal RNA (the anti-SD sequence). In some cases, internal SD-like motifs in the coding region generate "programmed" ribosomal pauses that are beneficial for protein folding or accurate targeting. On the other hand, such pauses can also reduce protein production, generating purifying selection against internal SD-like motifs. This selection should be stronger in GC-rich genomes that are more likely to harbor the G-rich SD motif. However, the nature and consequences of selection acting on internal SD-like motifs within genomes and across species remains unclear. We analyzed the frequency of SD-like hexamers in the coding regions of 284 prokaryotes (277 with known anti-SD sequences and 7 without a typical SD mechanism). After accounting for GC content, we found that internal SD-like hexamers are avoided in 230 species, including three without a typical SD mechanism. The degree of avoidance was higher in GC-rich genomes, mesophiles, and N-terminal regions of genes. In contrast, 54 species either showed no signature of avoidance or were enriched in internal SD-like motifs. C-terminal gene regions were relatively enriched in SD-like hexamers, particularly for genes in operons or those followed closely by downstream genes. Together, our results suggest that the frequency of internal SD-like hexamers is governed by multiple factors including GC content and genome organization, and further empirical work is necessary to understand the evolution and functional roles of these motifs.

Large-Effect Beneficial Synonymous Mutations Mediate Rapid and Parallel Adaptation in a Bacterium.

Contrary to previous understanding, recent evidence indicates that synonymous codon changes may sometimes face strong selection. However, it remains difficult to generalize the nature, strength, and mechanism(s) of such selection. Previously, we showed that synonymous variants of a key enzyme-coding gene (fae) of Methylobacterium extorquens AM1 decreased enzyme production and reduced fitness dramatically. We now show that during laboratory evolution, these variants rapidly regained fitness via parallel yet variant-specific, highly beneficial point mutations in the N-terminal region of fae These mutations (including four synonymous mutations) had weak but consistently positive impacts on transcript levels, enzyme production, or enzyme activity. However, none of the proposed mechanisms (including internal ribosome pause sites or mRNA structure) predicted the fitness impact of evolved or additional, engineered point mutations. This study shows that synonymous mutations can be fixed through strong positive selection, but the mechanism for their benefit varies depending on the local sequence context.